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The onset and evolution of magnetic fields in laboratory and astrophysical plasmas is determined 
by several mechanisms [T], including instabilities [2 El, dynamo effects HEl and ultra-high energy 
particle flows through gas, plasma and interstellar-media ElE]. These processes are relevant over 
a wide range of conditions, from cosmic ray acceleration and gamma ray bursts to nuclear fusion 
in stars. The disparate temporal and spatial scales where each operates can be reconciled by 
scaling parameters that enable to recreate astrophysical conditions in the laboratory. Here we 
unveil a new mechanism by which the flow of ultra-energetic particles can strongly magnetize the 
boundary between the plasma and the non-ionized gas to magnetic fields up to 10-100 Tesla (micro 
Tesla in astrophysical conditions). The physics is observed from the first time-resolved large scale 
magnetic field measurements obtained in a laser wakefield accelerator. Particle-in-cell simulations 
capturing the global plasma and field dynamics over the full plasma length confirm the experimental 
measurements. These results open new paths for the exploration and modelling of ultra high energy 
particle driven magnetic field generation in the laboratory. 
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Figure 1: a: Scheme of the experimental setup, showing the laser pump beam (red) focused into the gas jet and the probe 
beam (blue) split multiple times for simultaneous measurement of electron density and Faraday rotation in the plasma. The 
probe beam propagates perpendicularly to the pump laser path, integrating the polar magnetic component and the density in 
the plasma. b,c: Images of the plasma while the laser propagates at t = 1.9 ps. b: interference map showing the integrated 
phase, c: one of the two simultaneous images recording the polarization rotation onto the intensity (false colors). 

Strong cosmic magnetization requires ultra-high energy particle flows. These non-thermal particle streams can be 
originated [8] by statistical acceleration processes, such as shock 13 nm and Fermi acceleration m, and by direct 
particle acceleration mechanisms in strong wave fields, such as those found in pulsars m In addition to direct 
astronomical observations, the physical processes occurring when these mechanisms take place may also be explored 
in the laboratory. For instance, the conditions for the onset of statistical acceleration mechanisms through collisionless 
shocks in the laboratory have been investigated theoretically m- Direct cosmic acceleration can also be explored 
in the laboratory through laser driven plasma wakefields [TiHre]. Here we show that the non-thermal particle flows 
produced in a laser wakefield accelerator (LWFA) pTl [T5] can strongly magnetize the plasma and the plasma-neutral 
gas boundary. This observation is also the first time the strong magnetization occurring at the flow of energetic 
particles from ionised to non-ionised interstellar material can be reproduced in the laboratory. 
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A laser wakefield accelerator (LWFA) uses short and intense laser pulses to drive large amplitude plasma waves. 
In the LWFA scheme, the laser ponderomotive force excites ultra-relativistic, large amplitude plasma waves where 
electrons can be trapped and accelerated to high energies. In our experiment we used a 30 fs laser pulse focused to 
a transverse spot-size of 8 pm (FWHM) for a peak intensity of Iq = 3 x 10^^ W/cm^. The plasma is created in an 
Helium gas jet at near critical density, ua = 3.5 x 10^^ atoms/cm^. Since the laser power is above the critical power 
for self-focusing, it excites strongly non-linear plasma waves above the wave-breaking threshold. Wave-breaking leads 
to non-thermal particle flows with relativistic energies. 

Energetic particles produced during wavebreaking propagate mainly in the forward (laser) direction. However, a 
fraction of these particles also propagates radially to regions of undisturbed plasma, eventually reaching the plasma- 
neutral gas boundary. As they propagate through the ionised medium, return currents are set up to balance the hot 
electron flow, preventing effective magnetic field generation. However, as the hot electrons cross the plasma-neutral 
gas boundary, strong magnetic fields are induced by the resulting current imbalance. We demonstrate that this mech¬ 
anism can produce strong magnetic fields scaling with 32/^ hot (no[10^^cm“^])^^^ T in the laboratory. In astrophysical 

conditions, fields created by the same mechanism can reach amplitudes of (0.1 — 1) r^hot (^o[cm“^])^^^ pT, where r^hot 
is the fraction of hot electrons/energetic particles to the background plasma density no- 

Our experiment provides the first time resolved measurements of the magnetic field spatial distribution in a laser- 
plasma accelerator for the whole plasma volume and with a high temporal resolution. As the laser propagates through 
the gas target we observed strikingly complex magnetic structures, with several inversions of the field orientation at 
the plasma core. These observations are in excellent agreement with 3D one-to-one Particle-in-Cell (PIC) simulations 
in OSIRIS [20] capturing the global plasma dynamics and magnetic field evolution over the entire gas jet each 
taking several hundreds of thousands of CPU-hours (see methods section). 



Figure 2: Magnetization of the plasma at different times. The laser pulse enters the gas area at t = 0.3 ps from the left hand 
side of the box. t=0.7ps, the bell-shaped plot represents the longitudinal electron density profile in the plasma, highest density 
being rie = 7 x 10^^ cm“^. The grayed region represents the beam waist projection. 

The laser interaction with the gas jet is thoroughly scanned at high temporal and spatial resolutions for electron 
density (via phase recording) and magnetic field mapping. The plasma is probed by a single 30 fs pulse, synchronized 
to the pump pulse and doubled in frequency, as shown in Figj^a. During its propagation in the plasma, polarization 
rotation and absolute de-phasing are integrated and recorded on three separate CCD cameras. The three-dimensional 
density and magnetic maps are then reconstructed from the recorded information (see Methods). This pump-probe 
experiment permits us to “freeze” the plasma state to a high resolution snapshot lasting only for 30 fs and to follow its 
evolution. The validity of our diagnostic reposes on the assumption that B • Vne = 0 in the plasma cylinder, which 
is confirmed by simulations. 

Snapshots of the spatial distribution of the magnetic field at selected times are shown in Figure Each image 
represents the symmetrized radial map of (ngB) , (r, x), product between the polar component of the magnetic field 
and the local electron density, as reconstructed n*om the probe polarization maps. The laser pulse propagates from 
the left to the right. A polar magnetic field is observed in the trail of the laser pulse soon after its entrance in the gas 
jet (t = 0.7ps). Erom the density measurement we can infer a magnetic field magnitude reaching 100T. This field 
is positive in the plasma core (r < 50 pm) and changes in sign around the radial border of the plasma (r ^ 50 pm). 
As we will see from simulations below, the inversion of the field direction at the boundary is a distinctive feature 
indicating an electron current passing trough the plasma/gas boundary. These electrons (ultra-high energetic particles 
in astrophysical scenarios) are relativistic and, though accelerated from the wavebreaking, they are not trapped by 
the wakefield. 
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Figure 3: Composed images showing the evolution in time of the magnetic held (temporal resolution: 200 fs). The three panels 
correspond to a, b, c markers in Fig[^ Dotted lines indicate the laser arrival. 


Figure shows another striking feature, consisting in a strong magnetization of the plasma core. A large magnetic 
held is indeed expected close to the laser propagation axis, in the wakeheld region. This held is created by the very 
high longitudinal wakeheld current, which is not screened by the plasma. At a distance corresponding to the limit of 
the laser beam waist (i.e. a few pm’s), a weak magnetic held survives (the bow wave magnetic held [21]). Farther away 
from the axis, plasma return currents typically screen hot electrons that expand radially from the wakeheld (note 
that the main velocity component of these electrons is along x), thus no signihcant held is expected to survive up to 
the boundary whereabouts. The experimental results suggest that in the plasma core the screening of hot electron 
current coming from the wakeheld is not as efficient as expected. 

The magnetic held structure remains alike until t = 1.1 ps: starting from t = 1.3 ps, an island of reversed held starts 
building up close the gas jet density peak {x = 530pm), relatively far from the laser waist. At t = 1.5ps Q, this 
island has expanded longitudinally and radially, and a second one appears closer to the laser axis {x = 580 pi^, at a 
larger radius. Simulations suggest that these islands are caused by the hlamentation of the wings of the laser pulse, 
i.e. the laser energy surrounding the central spot. 

Experiments showed several inversions of the magnetic held orientation as illustrated in Fig. Several shots at 
different delays are used to compose a picture of the local evolution in time of the magnetic held. The temporal 
resolution in this hgure is of 200 fs. Upon its entrance in the gas (FigJ^) the laser propagation is accompanied by 
the formation of a positive poloidal magnetic held in the plasma core, which changes sign in the plasma boundary. 
The inner held component is consistent, in sign, with a negative current propagating with the laser pulse and remains 
stable for approximatively 1.2 ps. Panel b of hgure shows the evolution in the region where the hrst island appears, 
X = 480 pm. At t ^ 1.1 ps the held sign is reversed in the plasma core, corresponding to the formation of the hrst 
island. This inversion lasts up to t ^ 1.5 ps. When the laser exits from the gas prohle (FigJ^) a poloidal held is 
reformed, comparable to panel (a). At this time, however, the plasma core magnetisation is opposite, in sign, to what 
was observed at the entrance. 

Experimental hndings are conhrmed by 3D PIC simulations run with parameters closely matching experimental 
laser and plasma conditions (see methods). As it enters the gas, the laser ionizes the gas up to a radius of 100pm 
away from the axis and excites weakly non-linear plasma waves. Relativistic pulse self-focusing enhances the wakeheld 
amplitude beyond the wavebreaking threshold after 400 pm of propagation. When wavebreaking occurs, a fraction of 
the resulting hot electrons expands radially through the plasma. 

The laser driven plasma waves lead to complex longitudinal electron current structures in the ionised volume. 
Eigurej^I shows the current structures at the end of the simulation, at t = 2.67 ps, where it is possible to distinguish 
between backward electron currents (blue) and forward electron currentsf red) . As electrons cross the plasma-neutral 
gas boundary an inner return current is set up, located at the boundary itself, at a variable radius in the range 
r < 60 — 100 pm. These currents, hgure [^I, are at the origin of the large scale poloidal magnetic helds observed in 
hgure [^II. Simulations also show that a small fraction of hot electrons at the plasma entrance {x < 500 pm) how away 
from the laser in the backward direction. These are indicated by the red structures involving a blue bulk (r > 60 pm) 
in Eig. I^Ia. Eor x > 500 pm most of the hot electrons are accelerated in the forward direction leading to the outer 
forward (blue) current structures in Eig. §I b, c. 
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Return currents electrons propagate forward at the leading edge of the gas profile (Fig. la) and backward for 
the remainder of the gas jet length (Fig. [^Ib-c). Because hot electrons flow in opposite directions for x < 500 pm 
(Fig. Ila) and for x > 500 pm (Fig. |^IIc), the magnetic field sign changes in these regions. The transition at 
X 500 pm is shown in Fig. Ilb, which illustrates the poloidal magnetic field corresponding to the longitudinal 
currents shown in Fig. L This transition to a more rich current structure around x 500 pm is confirmed by 
experimental observation, e.g. Fig. and Fig. lit- 

The large scale magnetic field that surrounds the plasma around \y\ 60 — 100 pm reaches amplitudes that are in 

fair agreement with experimental measurements {n^Bz — 10^^ T/cm^). 



Figure 4: Summarized results of a full 3D PIC simulation of the experiment (see Methods). The final condition of currents (I) 
and azimuthal magnetic field (II) in the plasma is shown after the laser pulse has propagated through the gas; three 2D cuts 
are shown for improved readability. 

An estimate of the poloidal magnetic field amplitude can then be determined considering the return currents at 
the plasma-gas interface. Using Ampere’s law to estimate the amplitude of the resulting magnetic field in cylindrical 
symmetry gives /B • dl 27rrBo = (47r/c) /j • dS = po /j;s * dS, where r is the distance to the axis. For rela¬ 
tivistic hot electrons moving at c longitudinally j ecnoBz^ where e is the elementary charge and Bz is the unit 
vector pointing in the 2 ; direction. When r is larger than the plasma radius the return current density flux is 
^hot^nodTr^ [vp — {vp — A)^] ^ dTT^ceno 77hot^pA, where A c/uOp ^ r is the thickness of the plasma where the re¬ 
turn currents setup, hence Bq ^ AireriQ c77hot {c/ojp). The typical generated magnetic fields are on the order of 

-B0[r] ~ 327?hot\AoIi^ cm for r Vp. Considering 77 hot — 0-05 (taken from simulations) and no = 7 x 10 cm ^ 
we obtain Bq 80 T, which is consistent with the simulation results. In astrophysic al scenario s, the amplitude of the 
cosmic magnetic fields generated by this mechanism can be estimated as 0.32/7hot \/^o[cm“^] pT. 

The poloidal magnetic fields observed in our experiment also provide signatures for the occurrence of wavebreaking, 
which can be controlled by tuning laser and plasma parameters. The amplitude of the fields can be controlled by the 
plasma radius (or the typical size of the ionised cosmic region) and also by the laser intensity, which determines the 
fraction of hot electrons (or ultra high energy particles such as cosmic rays). The orientation of the field indicates the 
preferential direction of the hot electron flow. 


Methods 

High density micrometric target The target is a pulsed, high pressure, gas jet system which can drive a sub-millimetric 
gas nozzle to atomic densities in the range 10^^ — 10^^ atoms/cm^ [22] . 

In our experiment we used a transonic nozzle with an output diameter of 400 pm, producing an expanding flow with a Mach 
number M = 1.3. The radial atomic density profile at the distance of 200 pm from the nozzle level (laser propagation axis) is 
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fit by ua = no exp (r/ro)^^‘^^J, where no = 3.5 x 10^^ atoms/cm^ and ro = 170 pm is the jet radius. 

The peak density in the jet exponentially decreases along the vertical direction with the characteristic scale length of 263 pm. 
Probing of the plasma The plasma is probed by a linearly polarized Tp = 30 fs, A = 400 ns probe pulse synchronized to the 
main (pump) beam. The probe beam is split onto three separate optical setups, for simultaneous measurement of density and 
polar magnetic field. 

The radial density map is extracted from the integrated phase of the probe pulse, via a Nomarski interferometer on the first 
transmitted copy of the probe beam. The integrated phase is retrieved by 2D wavelet analysis of the interferogram. The phase 
map is then normalized by subtraction of a reference phasemap (accounting for aberrations in the laser transport or in the 
imaging system) and inverted by Hankel-Fourier implementation of the Abel transform. 

Magnetic held is retrieved from the polarization rotation induced on the probe beam (Faraday effect) as 


^rot — 


e 

2 me cric 


I 


rie B • ds 


( 1 ) 


where Uc = me^o (27rc/eA) is the electron plasma critical density at the wavelength A (see [23]). 

In order to obtain a 2D map of the polarization status in the probe beam prohle, the pulse is split onto two separate diagnostic 
lines, each equipped with an analyser and a high dynamic range CCD camera. A total of four images, two with pump and two 
without, are used for each snapshot of the magnetic held to eliminate ehects due to laser intensity huctuation and systematic 
optical deformation. The correct superposition of images from the two separate polarization measurement lines is ensured by 
spatial markers and automated numerical pattern recognition methods. 

Considering that the ratio map R{x,y) = Ipi/Ip 2 can be written as 


R (x, 

where ^p{i, 2 } are the analyzers angles, the local 
and Ap = 0 p 2 — 0pi from 

0i 


^ _ 1 - /3i sin^ ('k/2 + 8pi + iprot) 

1 — /32 sin^ ('^/2 + 0p2 + ^rot) 

polarization rotation (frot {x^ y) can be calculated dehning 0i — 7r/2 + ^pi if rot 


= tan 


a+{a'^ +h'^ 
5 + c 


(3) 


where coefficients are given by 


{ a — Rj32 cos (2Ap) — /3i 
b = —Ri32 sin (2Ap) 
c = RI32-I3i^2{1-R). 


(4) 


The 2D map of the quantity {ue is hnally obtained from the projected polarization rotation using the appropriate integral 
transform [24]: 


{neB^){r,x) = 


2mecnc d frot{y,x) 


dr 


nr 

J r 




dy. 


Due to the sensitivity of the algorithm, frot maps are antisymmetrized before back transformation. 


(5) 


Simulations Three-dimensional simulations were performed using the fully relativistic partice-in-cell code Osiris [191120] . which 
is routinely used to model laser wakeheld acceleration and astrophysical scenarios. The simulation window is 1272 x 300 x 
300 (pm)^, divided into 33000 x 1000 x 1000 cells with 2x1x1 particles per cell, giving a total of 6.6 x 10^° simulation particles 
and a total simulation time between 200000-400000 CPUh. Simulations considered an initial Helium gas jet density prohle given 
by no = 3.5 x 10^^ cm“^ exp [(| X — 700 pm| /340pm)^'^], reproducing experimental conditions. Ionisation was modelled using 
ADK tunnel ionisation rates. The longitudinal prohle of the laser electric held is symmetric and given by 10 — 15 t^ + 6t^, 

with T = \/2t/TFWHM, where tfwhm is the FWHM duration of the laser pulse. The transverse laser prohle is a ht to a transverse 
line out of the experimental laser intensity prohle using higher order Hermite-Gaussian modes. Each mode electric held prohle 
is dehned as: 


(^) Hm (^) • 


( 6 ) 

• (i^) + 2(x°+lp) - Cm,u(a:)] , 


where Eq is the peak electric held, ko = ujq/c, uq = 2.34 x 10^^ rads“^ is the central laser frequency corresponding to a central 
laser wavelength of 800 nm, Cp{x) = (m + n + 1) tan“^(x, Z^.) is the Gou^M^hase shift, = /coVFo/2 is the Rayleigh length, 
W{xy = VFo (1 -1- x^lZl\ and Hn the n-order Hermite polynomial. Tableshows the laser peak normalised vector potentials 
(ao) of the hve higher order Gaussian beams associated with the ht to the experimental laser prohle. In addition. Wo — 9.8 pm 
was also used. 
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(m,n) (0,0) (0,1) (0,2) (0,3) (0,4) 

Vo 2.27 -0.076 -0.05 6.3 x 10“'’ 0.025 


Table I: Laser intensities for the higher order Hermite-Gaussian modes used to fit the experimental laser profile. 
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